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I present a recent piece of work on semilcptonic B — » tt decays in which lattice results and kinematical and 
dispersive constraints are combined to obtain model-independent bounds on the relevant form factors and rates. 



1 Introduction 

The CLEO Collaboration has very recently pre- 
sented a measurement of -the branching ratios for 
B — > ixlv decays (£=e,/i)EJ: 

B (B° -> 7r+r^) = (1.8±0.4±0.3±0.2) x 1CT 4 , 

(1) 

where the errors are statistical, systematic and the 
estimated mo del- dependence introduced in deter- 
mining efficiencies. This measurement represents 
an excellent opportunity to determine the poorly 
know CKM matrix element \V u b\- This determina- 
tion requires calculating, over the full kinematical 
range, the non-perturbative matrix element: 

(w+(p')\V»\B°(p)) = (p + p'-qMl^y f+ 

+q« M ^f° , (2) 

where q—p—p', V^—u^^b, M is the mass of the 
B and m, that of the tt. Though this matrix ele- 
ment can be calculated using lattice QCD, cur- 
rent day lattice simulations, with lattice spac- 
ings on the order of 3 GeV -1 , cannot cover the 
full kinematical range. The problem is that the 
energies and momenta of the particles involved, 
whose orders of magnitude are set by the b quark 
mass (m& — 5 GeV), are large on the scale of 
the cutoff over much of phase space. To limit 
these energies in relativsitc lattice quark calcu- 
lations, one performs the simulation with heavy- 
quark mass values toq around that of the charm 
(m c ~ 1.5 GeV), where discretization errors re- 
main under control. Then one extrapolates the 
results up to nib by fitting heavy-quark scaling re- 
lations (HQSR) with power corrections to the lat- 
tice results. Another approach is to work with 
discrctizcd versions of effective theories such as 
Non-Rclativistic QCD (NRQCD) or Heavy-Quark 
Effective Theory (HQET) in which the mass of 
the heavy quark is factored out of the dynamics. 



All these approaches, however, are based on the 
heavy-quark expansion which is more limited for 
heavy — > light quark decays, such as the one that 
concerns us here, than it is for heavy — > heavy 
quark decays: it imposes no normalization condi- 
tion on the relevant form factors at the the zero 
recoil point q 2 = q^ax anc ^ oru y applies in a limited 
region around q^ax- Furthermore, momentum- 
dependent discretization errors restrict the mo- 
mentum of the initial and final state mesons to 
around 1 ~ 2 GeV. Thus, all these approaches are 
constrained to relatively small momentum trans- 
fers: one can only reconstruct the q 2 dependence 
of the relevant form factors in a limited region 
around q 2 nax and one is left with the problem of 
extrapolating these results to smaller q 2 . 

heavy — ► light quark decays are difficult 
in any theoretical approach. Indeed, they re- 
quire understanding the underlying QCD dynam- 
ics over a large range of momentum transfers from 
q 2 nax =26A GeV 2 for semileptonic B — > tt decays, 
where the tt is at rest in the frame of the B meson, 
to q 2 =0 where it recoils very strongly. 

2 B° — > tt + £~D and Dispersive Constraints 

The solution to the problem of the limited kine- 
matical reach of lattice simulations of heavy — > 
light quark decays presented here consists in sup- 
plementing lattice results for the relevant form fac- 
tors around q 2 nax with dispersive bound techniques 
to obtain improved, model-indgpendent bounds 
for the form factors for all q 2 . u For the case of 
B Q — > n + l~v decays, one can use the kinematical 
constraint, / + (0)=/°(0), to further constrain the 
bounds. 

2.1 Dispersive Bounds 

The subject of dispersive bounds in semileptonic 
decays has a long history going back to S. Okubo 
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et aLwho applied them to semileptonic K — > it de- 
caysS. C. Bourrely et al. first combined these tech- 
niques with QCD and applied them to semilep- 
tonic D — > K decaysQ. Very recently, C-G. Boyd 
et al. applied them to B — > irlv decays a. 

The starting point for B — > irlD decays is the 
polarization function 

W v {q)=i J d 4 x e iq - x {0\T(V»(x)V^(0)) |0) 

=W - </V ) nr(g 2 ) + n L tf) , (3) 

where Ht(l) corresponds to the propagation of a 
J p =l~ (0 + ) particle. The corresponding spectral 
functions, ImllT.L, are sums of positive contribu- 
tions coming from intermediate B* (J p =l~), Bit 
(J p =0 + and 1~), ... states and are thus upper 
bounds on the Bit contributions. Combining, for 
instance, the bound from Imlli with the disper- 
sion relation (Q 2 = — q 2 ) 

xl(Q 2 ) = ^(Q 2 n L (Q 2 )) 



dQ 2 
1 



dt 



tImIL L (t) 
(t + Q 2 f 



(4) 



one finds 



xUQ 2 ) > 



dtk(t,Q 2 )\f°(t)f 



where t±=(rriB ±TOtt) 2 and k(t,Q 2 ) is a kinemat- 
ical factor. Now, since xl{Q 2 ) can be calculated 
analytically in QCD for Q 2 far enough below the 
resonance region (i.e. — Q 2 <C m 2 ), Eq. ^ gives 
an upper bound on the weighted integral of the 
magnitude squared of the form factor f° along the 
Bit cut. To translate this bound into a bound on 
f° in the region of physical B — > n£p decaysJs a 
problem in complex analysis (please see Ref.13 for 
details). A similar constraint can be obtained from 
ILt for / + . There, however, one has to confront 
the additional difficulty that / + is not analytic 
below the Bit threshold because of the B2 pole. 

The beauty of the methods of Ref. □ is that 
they enable one to incorporate information about 
the form factors, such as their values at various 
kinematical points, to constrain the bounds. For 
the case at hand, however, these methods must be 
generalized in two non-trivial ways. In construct- 
ing these generalizations, one must keep in mind 
that the bounds: 1) form inseparable pairs; 2) do 
not indicate the probability that the form factor 
will take on any particular value within them. 



2.2 Imposing the Kinematical Constraint 

The first problem is that Eq. (|^) and the equiv- 
alent constraint for / + yield independent bounds 
on the form factors which do not satisfy the kine- 
matical constraint / + (0)=/°(0). The bounds on 
/ + require / + (0) to lie within an interval of val- 
ues J+ and those on /°, within an interval Iq. To- 
gether with these bounds, however, the kinemat- 
ical constraint requires / + (0)=/°(0) to lie some- 
where within 7+nIo. Thus, we seek bounds on the 
form factors which are consistent with this new 
constraint. 

A natural definition is to require these new 
bounds to be the envelope of the set of pairs of 
bounds obtained by allowing / + (0) and /°(0) to 
take all possible values within the interval i+ (~l 
lo- In Ref. □, it is shown how this envelope can 
be constructed efficiently and that the additional 
constraint can only improve the bounds on the 
form factors for all q 2 . Also, as a by product, one 
obtains a formalism which enables one to constrain 
bounds on a form factor with the knowledge that 
it must lie within an interval of values at one or 
more values of q 2 . 



(5) 2.3 Taking Errors into Account 



As they stand, the methods of Ref. □ can only 
accommodate exact values of the form factors at 
given kinematical points and contain no provisions 
for taking errors on these values into account. Of 
course, the results given by the lattice do carry 
error bars. More precisely, the lattice provides a 
probability distribution for the value of the form 
factors at various kinematical points. What must 
be done, then, is to translate this distribution into 
some sort of probability statement on the bounds. 
The conservative solution is to consider the prob- 
ability that complete pairs of bounds lie within 
a given finite interval at each value of q 2 . Then, 
using this new probability, one can define upper 
and lower p% bounds at each q 2 as the upper and 
lower boundaries of the interval that contains the 
central p% of this probability.^ These bounds in- 
dicate that there is at least a p% probability that 
the form factors lie within them at each q 2 . 



"The density of pairs of bounds increases toward the 
center of the distribution as long as the distribution of the 
lattice results does. 
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2-4 Lattice- Constrained Bounds 

To constrain the bounds on / + and fP, the lattice 
results of the UKQCD Collaboration^] are used, to 
which a large range of systematic errors is added to 
ensure that the bounds obtained are conservative. 
Because of these systematic errors, the probabil- 
ity distribution of the lattice results is not known. 
The simplifying and rather conservative assump- 
tion that the results are uncorrelated and gaussian 
distributed is made. The required probability is 
constructed by generating 4000 pairs of bounds 
from a Monte-Carlo on the distribution of the lat- 
tice results. The results for the bounds on the 
form factors are shown in Fig. [I]. The two form 
factors are plotted back-to-back to show the ef- 
fect of the kinematical constraint. Without this 
constraint, the bounds on /+ would be looser, es- 
pecially around q 2 —Q, where phase space is large. 
Since / + determines the rate, the kinematical con- 
straint and the bounds on /° are important. 

Also shown in Fig.Ij] is the light-cone sumrule 
(LCSR) result of Ref. El and the 3-point sumrule 
(SR3) result of Ref.13. The latter is parametrized 
by a poleJorm with / + (0) = 0.26 and m pole — 
5.25 GeVcl. The former has two components: for 
q 2 below 15 GeV 2 , the q 2 dependence of /+ is 
determined directly from the sumrule; for larger 
q 2 , pole dominance is assumed with a residue de- 
termined from the same correlator. While agree- 
ment of the LCSR result with the bounds is excel- 
lent, the SR3 result is quite strongly disfavored. 
The bounds are compared with the predictions of 
more authors as well as with direct fits of vario 
parametrizations to the lattice results in Ref. 
Though again certain predictions are strongly dis- 
favored, the lattice results and bounds will have to 
improve before a firm conclusion can be drawn as 
to the precise q 2 dependence of the form factors. 

The bounds on / + also enable one to con- 
strain the B*Bn coupling gs'Bn which determines 
the residue of the B* pole contribution to /+. 
The constraints obtained are poor because /+ is 
weakly bounded at large q 2 , as can be seen in 
Fig. yj. Fitting the lattice results for f° and / + 
to a parametrization which assumes B* pole dom- 
inance for / + and which is consistent with HQS 
and the kinematical constraint gives the more pre- 
cise result gB'+B°TT+ — 28 ± 4.|| However, because 




Figure 1: f°(\q 2 \) and f + (q 2 ) vensus q 2 . The data points 
are the lattice results of UKQCD □ with added systematic 
errors. The pairs of fine curves are, from the outermost to 
the innermost, the 95%, 70% and 30% boiinids. The upper 
shaded curve is the LCSR result of Ref. fc and the lower 
one, the SR3 results of Ref. Q. 



this result is model-dependent, it should be taken 
with care. 

2.5 Bounds on the rate and \V u b\ 

As was done for the form factors, one can de- 
fine the probability of finding a complete pair of 
bounds on the rate in a given interval and from 
that probability determine confidence level (CL) 
intervals for the rate. The resulting bounds are 
summarized in Table [l]. They were obtained by 
appropriately integrating the 4000 bounds gener- 
ated for /+ (q 2 ) , taking the skewness of the result- 
ing "distribution" of bounds on the rate into ac- 
count. Also given in Table |l| are bounds on / + (0) 
as well as the preditions of other authors for these 
quantities. Although the bounds on f + (q 2 ) disfa- 
vor some of these predictions rather strongly, the 
bounds on the_rate only exclude the quark model 
result of Ref. E3 with high certainty. 

The CL bounds of Table [l] can also be used, in 
conjunction with the branching ratio measurement 
of CLEO given in Eq. (|]), to determine \ V u b\. One 
finds 



6 The result of this fit is entirely compatible with our 
bounds on /+ and f°. 



\V ub \10 4 y/T B o /1.56ps = (34-f49)±8±6 , (6) 

where the range given in parentheses is that ob- 
tained from the 30% CL bounds on the rate 
and represents the most probable range of val- 
ues for \V u b\- The first set of errors is obtained 
from the 70% CL bounds and the second is ob- 
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Table 1: Bounds on 
/+(0): top block. 
(LCSR), 2pt (SR2) 



rate in units of \V u i,\ ps and on 
Quark model (QM) and light-cone 
and 3pt (SR3) sumrule predictions: 
bottom block. 



1 (is — > ir^Z v\ 


f (0) 


details 


2.4 — » 28 


—0.26 — * 0.92 


95% CL 


2.8 -> 24 


-0.18 -> 0.85 


90% CL 


3.6 — > 17 


0.00 — ► 0.68 


70% CL 


4.4 — > 13 


0.10 — > 0.57 


50% CL 


4.8 — > 10 


0.18 — > 0.49 


30% C 


L 


7.4 ± 1.6 


0.33 ± 0.06 


Qlvd 




2.1 


0.09 


QM L 




9.6 




QML 


i 




0.26 


QM L 






0.24-0.29 


QML 




9.6-15.2 


0.29 - 0.46 


QM L 




7±2 


0.20-0.29 


QML 




14.5 ±5.9 


0.4 ±0.1 


SR2tL 


i 


4.5-9.0 


0.27 ±0.05 


SR2 +1 


> S 


3.60 ±0.65 


0.23 ±0.02 


SR3E. 


8 


5.1 ±1.1 


0.26 ±0.02 


SR3t 




8.1 


0.24-0.29 


LCSR 




7± 1 


0.21 - 0.27 


latJ 




9±6 


0.10-0.49 


LAT I 


9 


8±4 


0.23-0.43 


LATE 


() 



tained by combining all experimental uncertain- 
ties in quadrature and applying them to the av- 
erage value of \V u b\ given by the 30% CL results. 
This determination of \V u b\ has a theoretical er- 
ror of approximately 37%. Though non-negligible, 
this error is quite reasonable given that the bounds 
on the rate arc completely modcl-indcpcndcnt and 
are obtained from lattice data which lie in a lim- 
ited kinematical domain and include a conserva- 
tive range of systematic errpis. For comparison, 
the lattice results of ELC E3 and APE l°l were 
obtained from the determination of / + at a sin- 
gle q 2 ~ 18 — 20 GeV 2 supplemented with the 
assumption of £?*-pole dominance while those of 
UKQCD □ rely on constrained model fits to the 
same values of / + and f° (without the added sys- 
tematics errors) used to obtain the bounds pre- 
sented here. 

3 Conclusion and Outlook 

I have presented a new formalism by which lattice 
results for semileptonic B — > n decays, limited to 
a narrow kinematical range, are combined with 



dispersive and kinematical constraints to obtain 
model-independent bounds on the relevant form 
factors and rates over the full kinematical domain. 
I have compared these bounds, which have a well 
defined statistical meaning, to the predictions of 
other authors. 

Though the bounds will benefit from forth- 
coming, improved lattice results, they would ben- 
efit most from an increase in the range of q 2 . 

Finally, the techniques presented here are in 
principle applicable to limited results obtained by 
non-lattice means as well as to other processes 
such as B — * plv and B — > A"*7 decays. 
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